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ABSTRACT 

We utilize the extensive datasets available for the Perseus molecular cloud 
to analyze the relationship between the kinematics of small-scale dense cores 
and the larger structures in which they are embedded. The kinematic measures 
presented here can be used in conjunction with those discussed in our previous 
work as strong observational constraints that numerical simulations (or analytic 
models) of star formation should match. We find that dense cores have small 
motions with respect to the 13 CO gas, about one third of the 13 CO velocity 
dispersion along the same line of sight. Within each extinction region, the core- 
to-core velocity dispersion is about half of the total ( 13 CO) velocity dispersion 
seen in the region. Large-scale velocity gradients account for roughly half of 
the total velocity dispersion in each region, similar to what is predicted from 
large-scale turbulent modes following a power spectrum of P(k) oc k~ 4 . 



INTRODUCTION 



Star formation research is entering an era where the properties of dense, star-forming 
cores and their environs can be described in a statistical manner across many molecu- 
lar cloud environments. Existing multi -cloud surveys inc lude the COMPLETE Survey 



(IRidge et al.l 120061 ) . Spitzer c2d Survey (jEvans et al.ll2003l ). and Spitzer Gould Belt Sur- 
veys, while surveys that are in progress include the JCMT and Herschel Gould Belt sur- 
veys (j Ward- Thompson et al.l 120071 ; I Andre fc Saracenol 120051 ) . These surveys will character- 
ize dense core properties with sufficient detail to allow for the measurement of both global 
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properties as well as the variation due to molecular cloud environment. Numerical simula- 
tions have also begun to make explicit predictions about a variety of observable dense core 
properties, going beyond global statistics such as the mass function. 

One nearby molecular cloud that has been particularly well-studied already is P erseus. 
The Perseus molecular cloud is located relatively close-by at a distance of 250 pc (Cernis 
19931 ) and is an intermediate-mass star-forming region with protostars and dense cores found 
primarily within small clustered environments (e.g., NGC 1333 and IC 348). To study the 
dense cores in an unbiassed manner, dust column density observations have been undertaken 
in the (sub) millimetre continuum (Hatchell et al. 2005; Kirk, Johnstone, & Di Francesco 
2006, hereafter KJD06; Enoch et al. 2006) and the large-scale structure of the cloud has been 
determined through extinction inferred via reddening of stars in the 2MASS dataset (KJD06; 
Ridge et al. 2006). The young protostellar content of the cloud has been measured using 
Spitzer data (IRebull et al.ll2007l ; IjQrgensen et al.1120071 ; lHatchell et al.ll2007al ). Molecular line 
data also exist on multiple scales within the clo ud. The kinematics of the dense cores have 
been traced by poi nted observations of N^ H 4 " (IKirk. Johnstone fc Tafallal 120071 . hereafter 
KJT07) and NH -j (IRosolowskv et al.l 120081). and th eir immediate less-dense surroundings 
in C 18 (IKJT071 ) and CCS (IRosolowskv et al.l 120081 ) . Note that carbon-bearing molecules 
tend to freeze out quickly onto dust grains in the centres of cold, dense cores, while the 
absence of CO m olecules leads to a n enhanced production of nitrogen-bearing molecules 
(see, for example. iTafalla et al.ll2002l ). Hence in dense cores, nitrogen-based molecules tend 
to trace the dense parts of the cores, while carbon-bearing molecules tend to trace the 
lower density envelopes of the cores. The less-dense gas in Perseus has been mapped in its 
entire ty at lower resolution in 12 CO(1-0) and 13 CO(1-0) (IRidge et al.l 120061 ; IPineda et al. 



20081 ). with outflo ws identified acr oss the entire cloud in an unbiassed manner using the 
12 C O(1-0) dataset ( Arce et a and outflows around the dense cores mapped in 12 CO (3- 

2) jHatchell et alil2007bh . 



Since a wealth of data exist for the Perseus molecular cloud, it is possible to move 
beyond characterizing properties of structures in the cloud based on a single tracer, and 
advance to a study of the kinematic relationship between the structures on different scales 
within the cloud. 

Several studies have alre ady been carried out wh ich compare the column density features 
on various scales in Perseus. IjQrgensen et al.l (120071 ) related the young embedded protostars 
to the SCUBA dense cores that they inhabit and showed that the protostars have little, if 
any, motion within the dense core they form out of. In turn, the SCUBA dense cores are 
found only within high extinction zones within the mo lecular c loud, suggesting an extinction 
threshold of Ay ~ 5 mag for dense core formation (jKJD06f ). Such a threshold is consis- 
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tent with predictions from a magnetically-dominated formation scenario, wherein dense core 
formation is only possible in the highest column density regions, wher e the ionizatio n frac- 
tion is low enough to allow for a fast ambipolar diffusion timescale ( IMcKed Il989l ). The 
dense cores also tend to be offset from the peaks of the larger exti nction re gions that they 
inhabit in a pattern suggestive of a triggered formation scenario (1KJD06I ). A nearby BO 
star, 4 Persei, previously sugg ested as a trigger for star formation in the nearby L1451 
region (fWalawender et al.l 120041 ). is located in a position to be consistent with this scenario 
jKJDod ). 



A multi-scale analysis of kinematics comparable to the column density studies has not 
yet been performed, although the previous single-scale analyses reveal information which 
already places constrain ts on th e star formation scenario. The pointed N 2 H + (l-0) and 
C 18 0(2-l) observations (IKJT07I ) of the dense cores showed that the dense cores in Perseus 
are intimately tied to their surrounding envelope of material. The velocity dispersion of each 
core's N2H + gas tends to be thermally-dominated, while the velocity dispersion of each core's 
surrounding C 18 gas tends to be much more non-thermal. Despite this, the difference in 
centroid velocities between the N 2 H + and C 18 gas is typically less than the sound speed, 
and much less than the velocity dispersion of the C 18 Q gas. Similar resu lts were found from 



2008). In Perseus, 



the pointed NH 3 (1,1) and (2,2) and CCS(2-1) survey (IRosolowskv et al 
N2H + and NH3 have been shown to trace each other very well (1 Johnstone et al.ll2010l ) 



In a survey of isolated dense cores in a variety of star-forming regions, I Walsh et al. 



( 120041 ) found similarly small veloci ty differences and a rgued that this is evidence against 
the competitive accretion scenario. lAyliffe et al.l ( 120071 ) analyzed a simulation where cores 
under go compe titive accretion and produced mock observations to match the I Walsh et al. 
( 120041 ) survey. lAyliffe et al.l (120071 ) found they could match the small velocity differences 
at later times in the simulation, however, the dense core line width at later times greatly 
exceeded that found in observations. This illustrates that 'observations' of simulations are 
ne cessary for comparisons, as they can produce unexpected results. The Perseus observations 
of IKJT07I extend the result of small centroid velocity differences to a more clustered star 
formation en yironment, where competitive accretion is predicted to be more effective (see, 
for example, lAyliffe et al.l 120071 ). 



Using estimates of the mass and size of the dense cor es from the SCUBA observations 
and the internal velocity dispersion from the N2H + survey, IKJT07I furthermore showed that 
each dense core appears to be in approximate "virial equilibrium", so long as the surface 
pressure on the cores (from the overlying mass of the cloud) is considered. 

Simulators have started to 'observe' similar kinematic measures in their simulations, 
such as the internal velocity dispersion of cores (Klessen et al. 2005; Ayliffe et al. 2007; Offner 
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et al. 2008; Kirk, Johnstone, & Basu 2009, hereafter KJB09) and the core-to-envelope velocity 
difference (Ayliffe et al. 2007; KJB09). These measurements allow for a direct comparison 
between the observational kinematic surveys and the simulations. In this paper, we use 
the multi-scale kinematic data available in the Perseus molecular cloud to provide further 
observational measures by which to test future models of star formation. In the interest of 
facilitating the comparison between observations and simulations, we focus on physically- 
motivated quantities that are simple to derive. It is hoped that these simple observational 
measures can be related to simulations without many additional assumptions required for 
the simulations. 

In Section 2, we summarize the observational data used in our analysis, while in Sections 
3-6, we present results based on small- (Section 3) to large- (Sections 5 and 6) scale quantities. 
Interpretation of the results is discussed in Section 7, and we conclude in Section 8. A glossary 
of the kinematic measures used in this paper is provided in Appendix A. 



OBSERVATIONS 



Our analysis here is based primarily on two datasets - a pointed survey of the kinemat- 
ics of the dense cores (N2H + ) and their surrounding material (C 18 0) in Perse us k.IT07h. 
and a survey of the less -dense gas ( 13 CO) spanning the entire Perseus cloud (IRidge et al. 



20061 ; iPineda et al.ll2008l ). Additionally, we make use of the large structures (termed 'extinc- 
tion regions' in this paper) identified in an exti nction m ap of Perseus based on analysis of 
the 2MASS dataset ('extinction super cores' in IKJD06I ). Here, we summarize the relevant 
information for each of these datasets. 



2.1. Pointed N 2 H+ and C ls O Observations 



Using the IRAM 30 m telescope, IKJT07I made single pointing observations in frequency- 
switching mode of N 2 H + (l-0) and C 18 0(2-l) in 157 candidate dense cores in Perseus. The 
beamsize was 25" for N 2 H + and 11" for C 18 0. The pos itions observed were selected primarily 
based on the SCUBA submillimetre analysis of lKJDQq . with additional targets selected based 
on a visual inspection of the Palomar plates and the location of the m axima of large-scale 
structures seen in the 2MASS-based extinction map of Perseus (IKJD06I ). N 2 H + emission was 
detected in the majority of pointings (97), with the highest fraction of detections belonging 
to the SCUBA dense core list (84%) and lower fractions (42% and 14% respectively) for the 
others. C 18 emission was detected in all but one of the pointings. Linewidths, centroid 
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velocities, and integrated intensities were measured for all of the spectra with detected signal. 
N 2 H + (l-0) has seven hyperfine components in its spectrum and these were simultaneously 
fit using CLASS'S hyperfine fitting program. The vast majority of the N 2 H + spectra required 
a single velocity fit, and of the few cases where a second velocity component was required, 
the two components were usually discernible in the spectrum. The C 18 spectra were fit 
with a single Gaussian where appropriate, but more frequently required a two-component 
Gaussian fit. Table [I] summarizes the relevant results used in this analysis: the difference 
in centroid velocities, vn 2 h+ — vciso, an d the linewidth (one sig ma of a Gaussian fit) of the 



C 18 0, &ci80, for every N 2 H + detection^]. See Tables 3 and 4 in IkJT07| for the full spectral 



fit parameters for N 2 H + and C 18 respectively. Note that when N 2 H + spectra are plotted 
in this paper, the single Gaussian profile derived from the model fit is shown for simplicity. 

The dense cores were further divided int o protostellar an d starl ess, based on associations 
with the Perseus embedded YSO catalog of IJorgensen et al.l (120071 ). This YSO catalog uses 
two criteria of which protostars must satisfy at least one. The first criterion is the spatial 
coincidence of a Spitzer source and a SCUBA core, and the second criterion is colours that 
satisfy [3.6] - [4.5] > 1 and [8.0] - [24] > 4.5 in the Spitzer ba nds, i.e., colours that have been 
previously shown to select YSOs (see IJorgensen et al.l 120071 for more details). 



2.2. "CO Map 

The 13 CO(1-0) data were obtained by the COMPLETE teaiifl using the SEQUOIA 
focal plane array at the FCRAO telescope (IRidge et al.l 120061 ) . The region mapped covers 



~20 deg 2 , with a 46" angular resolution and 0.066 km s _1 velocity channels. The data 
were baseline-subtracted and sampled onto a 23" grid; the average noise is 0.1 K (in T* A 
sc ale). Detailed comp ariso ns between the molecul ar line data and extinction are presented 



in 



Pineda et al.l (|2008|) and iGoodman et al.l fl2009l ). 



In our comparison of the N 2 H + , C 18 0, and 13 CO spectra at the position of each core 
(Section 3), we fit a single Gaussian to each relevant 13 CO spectrum. The difference in 
N 2 H + and 13 CO velocity centroids and 13 CO velocity dispersions are given in Table [TJ All 
fits were visually examined; poor fits (i.e., where the spectra are not well- described by a 
single Gaussian) are also noted in Table 1. The estimation of the kinematic properties of 



1 A summary of all of the kinematic measures and their associated symbols used in this paper can be 
found in Table 5. 

2 AU of the Coordinated Molecular Probe Line Extinction Thermal Emission (COMPLETE) Survey data 
are publicly available at |http://www.cfa . harvard . edu/ COMPLETE/ 
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the larger extinction regions using the 13 CO data are discussed in Section 4. 



2.3. 2MASS Extinction 



Extinction maps are now readily available for the entire sky (see iRowles Froebrich 



20091 ). allowing large-scale structures to be easily identified in all molecular clouds. A map 
of the extinction throughout the Perseus cloud has been produced based on reddening of 
stars in the 2MAS S catalog (Ridge et al. 2 006; KJD06). This ma p was c reat ed using the 



NICE R technique (ILombardi & Alvesl l200ll ) , and was presented in IKJD06I and iRidge et al. 



(120061 ). The extinction map provides a useful dataset for defining the large-scale structures 
in Perseus. While this is not the only way to classify the large-scale column density struc- 
tures, it has the advantage of being simpler than disentangling structures in a spectral cube. 
In datasets with clustered emission features, which are particularly prevalent in 13 CO maps 
of molecular clouds (and generally more common in 3D datasets), the emergent properties 
of the structures identified are particularly sensi tive to the precise structure-identifying al- 
gorithm and parameters used (see, for example . IPineda et al.ll2009l ). For our analysis here, 
we use the large-scale structures identified in KJD06I using the 2D version of the clfind 
( [Williams, de Geus. fc Blitzl Il994f ) algorithm with a threshold of 3 mag and a stepsize of 
1 mag on the extinction map smooth ed to 5' resolution. T he releva nt properties of these 
extinction regions are summarized in iKJDpd (Table 3) and IkJT07I (Table 5 ). As a test 



we als o ran our analysis using the independently-defined large-scale regions in IPineda et al. 



(120081 ) based on CO centroid velocities and dust temperatures derived from IRAS data. Our 
conclusions are unchanged using these alternate regions, suggesting that our results are not 
strongly dependent on the definition of the regions used. We prefer to define the regions 
based on extinction, however, as this definition should be easy to apply to other observations 
and simulations. 

Figure [1] summarizes the data analyzed in this paper. The greyscale image shows the 
integrated intensity of the 13 CO data, while the red circles show the N 2 H + pointings. The 
total column density, as measured by the extinction observed in the 2MASS dataset is 
overlaid as black contours and the extinction regions are shown as coloured contours. 



RELATIVE MOTIONS AT THE DENSE CORE POSITIONS 



KJT07I measured the difference in centroid velocity of N 2 H + and C 18 (vn 2 h+ — Vciso), 
reflecting the motion between the dense cores (N 2 H + ) and their surrounding envelopes (C 18 0) 
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along the line of sight. This relative motion was found to be smaller than the envelope 
velocity dispersion (crciso) an d for 90% of the cores, less than the sound speed, assuming 
a temperature of 15 Ko- The small difference in centroid velocities implies that the dense 
core is dynamically coupled to its surrounding envelope, and appears to argue against the 
competitive accretion scenario (KJT07; Walsh et al 2004, but see also Ayliffe et al 2007). 
Since 13 CO(1-0) is sensitive to lower density material than C 18 0(2-l), we can additionally 
compare the N 2 H + and C 18 centroid velocities to that of the 13 CO to determine whether 
the tight correlation in velocities extends to the less dense material. (Note that the beamsize 
of the 13 CO observations is several times larger than that of the other tracers, so a larger 
physical volume of material is being traced, both in extent across the plane of the sky and 
along the LOS.) 

Figure |5] shows the distribution of centroid velocity differences as a function of RA (left 
panel) and the overall distribution of values (right panel). The filled symbols (left) and solid 
lines (right) denote locations where the 13 CO Gaussian fits were deemed 'good' (roughly 
87%), while the empty symbols and dotted lines denote locations where only the N 2 H + and 
C 18 fits were good. Note that the 13 CO map did not extend to L1455, hence no 'good' fits 
in that region were possible. All of these data are also provided in Table [TJ The greyscale 
image in the figure background shows the position-velocity (PV) diagram for the 13 CO gas, 
with the contours indicating the FWHM values. As can be seen from the figure, the spread in 
velocity difference is much smaller between N 2 H + and C 18 than between N 2 H + and 13 CO, 
and both of these are smaller than the range spanned by 13 CO emission at that RA. The 
standard deviation of the centroid velocity differences are 0.18 km s -1 for vn 2 h+ — vciso a ^ 
all locations, 0.19 km s _1 for VN2H+ — vciso where good 13 CO fits exist, and 0.33 km s" 1 for 
VN2H+ — vi3co , while the errors in the centroid velocity differences measured are roughly 
0.015 km s" 1 for N 2 H + and either CO species. A two-sided KS testo shows that the likelihood 
that both distributions are drawn from the same parent sample is only 2%. The sound speed 
is 0.23 km s" 1 for the mean gas, assuming a temperature of 15 K and a mean molecular 
weight of /i = 2.33. The largest velocity differences are found in the NGC 1333 region. 

A second way to measure the relative motion between the species is to compare the 
difference in velocity to the velocity dispersion of the less dense species. This normalized 



3 In terms of the full three-dimensional picture, the full (3D) velocity difference between the core and 
envelope would be expected to be larger than we observe, since only the LOS components of the velocities 
are measurable. Similarly, the envelope velocity dispersion and sound speed discussed also only reflect LOS 
components of the full 3D motion; all are equally affected by projection. 

4 A two-sided KS test allow s one to meas ure the probability of two distributions being drawn from the 
same parent distribution, (e.g.. lLuptonlll983l ). 
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velocity difference is written as: 

VN2H+ — VCO \ 

^norva 

&CO 

This quantity indicates whether the dense core moves within the typical range of motions 
present in the lower density en vironme nt, in this case, characterized by the gas along the 



same line of sight as the core. IKJB09I showed that the normalized velocity difference can 
provide strong discrimination between numerical simulations with varying input levels of 
turbulence and magnetic field strengths. The normalized velocity difference was also one 



of the most discrepant measure ments between many of the simulations analyzed in IKJB09 



and the observations of IKJT07I . Figure [3] shows the normalized velocity difference between 
the three species, as a function of RA (left panel) as well as the full distribution (right 
panel). As with the centroid velocity difference, we find the normalized velocity difference 
also tends to be quite small - nearly always less than one. The standard deviations of the 
normalized velocity differences are 0.20 (N 2 H+ - C 18 at all positions), 0.22 (N 2 H+ - C 18 
where good 13 CO fits exist), and 0.39 (N 2 H+ - 13 CO), while the errors are 0.014 and 0.019 
for N 2 H + to C 18 (all positions) and N 2 H + to 13 CO respectively. A two-sided KS test shows 
that the likelihood that both distributions are drawn from the same parent sample is 0.2%. 
This indicates that the motions of the dense and less-dense gas are correlated (i.e., closer 
than expected from a random distribution), but less so at larger scales. This is perhaps not 
surprising, since the 13 CO observations trace not only a larger scale across the plane of the 
sky (due to the larger beamsize), but also are expected to trace a much larger scale along the 
LOS, and likely include emission from additional lower-density structures not associated with 
the dense core, due to 13 CO's sensitivity to lower density material than C 18 0. Comparison 
with the results in Section 4.2 suggests that the agreement between N 2 H + and 13 CO centroid 
velocities might become slightly stronger if the 13 CO beamsize was small enough to match 
the N 2 H + beam. We expect, however, that the bulk of the difference is due to the larger scale 
of material traced by the 13 CO along the LOS, which would be unchanged with a smaller 
beamsize. 

A slight skew is apparent in both N 2 H + - 13 CO distributions, but does not appear to be 
physical. The mean velocity difference between N 2 H + and 13 CO is 0.03 km s _1 , less than half 
the size of the 13 CO spectral resolution, and hence is likely due to stochastic effects. As a 
second check, we examined the velocity differences within each extinction region and found 
that there was no systematic offset; offsets of opposite directions were found in different 
regions. 
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4. DENSE CORES RELATIVE TO THE EXTINCTION REGIONS 

4.1. Calculating Extinction Region Properties 
Since a full map exists for the 13 CO data, it is also possible to quan tify the motions of 



the dense cores within the larger extinction regions identified in iKJDOa In order to more 
easily identify large-scale features in the regions, we first convolve and re-grid the 13 CO 
map to match the 2MASS extinction map angular resolution (5') and grid. We create the 
mean 13 CO spectrum of the extinction region by s umming all 13 CO spectra that lie within 



the contours of the extinction region boundaries of IK JD06I . We then calculate the centroid 
velocity, v i3co,reg, and velocity dispersion, a r eg,Gauss, of the spectrum using a single Gaussian 
fit. Note that these spectra are little affected by the prior map convolution to 5'- in most 
cases, there is less than 3% variation between the Gaussian fit parameters for the total 
spectra of the conv olved an d un-convolved maps. Several of the spectra are not well-fit by a 



single Gaussian; in IKJB09I . we found that the full width of emission was better characterized 
by the full width at quarter maximum (FWQM). To determine the reliability of the Gaussian 
fits to the spectra, we also m easure the FWQM for the 13 CO region spectra, following the 



procedure outlined in IKJB09I . and calculate <J re g,FWQM, the sigma (1/e width) that would 
be measured if the spectrum were a perfect Gaussian with the measured FWQM. Table [2] 
provides the values fit to the cumulative spectra for each region using both methods, while 
Figure H] shows both results for all of the regions. Three of the eleven extinction regions have 
significantly incomplete 13 CO coverage. Regions that are not well-covered may not provide a 
good representation of the overall region kinematics, regardless of the method used to fit the 
total spectrum. We set a minimum value of 80% coverage in the 13 CO observations for all of 
our subsequent analysis. Figure @] shows that when a cutoff of 80% is applied, both measures 
of the total velocity dispersion (cr regi G auss or (J T eg,FWQM) produce similar results. For the 
remainder of the analysis and for the normalized velocity difference between the N 2 H + cores 
and their associated extinction regions provided in Table [U only the Gaussian fit results 
are used. Note that the tendency for the a T e g ,FWQM measures to be slightly lower than 
the corresponding a r eg,Gauss value appears to be primarily caused by our limited resolution. 
While the FWQM is a good measure for identifying fainter, broad components of spectral 
lines, when significant line wings are not present, noise and the small number of spectral 
channels will tend to make the FWQM underestimate the true width, as this measure is not 
interpolated across spectral channels (unlike the Gaussian fit). 
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4.2. Analysis 

Figure |5] shows the distribution of the centroid velocity differences between the N2H + and 
the 13 CO region, VN2H+ — ^xzco,reg versus RA (left panel), in addition to the full distribution 
of velocity differences (right panel). The motions between the core and region again tend to 
be small, although larger than the difference using the 13 CO velocity at the location of the 
dense core. The standard deviation of the distribution is 0.48 km s _1 (c.f. 0.33 km s _1 for 
13 CO gas at the dense core location); the error in the velocity difference measured is roughly 
0.011 km s -1 . Using a two-sided KS test shows that the probability of the distribution of 
velocity differences between the core and the 13 CO gas at the core's location versus over the 
entire region is only 0.5%. Figure E] shows the distribution of normalized velocity difference 
between the cores and regions, Cnorm,reg (calculated using eqn 1), versus RA (left panel) as 
well at the total distribution (right panel). The standard deviation of Cnorm,reg is 0.5 (c.f. 
0.39 for 13 CO gas at the dense core location) ; the error in the normalized velocity difference 
measured is 0.013. The two-sided KS test gives a probability that this normalized velocity 
difference and the one measured for the 13 CO at the core's location are drawn from the same 
parent distribution of 7%. We also examined the relationship between the velocity difference 
and the separation from the centre of the region, but found no discernable trend. 



5. CORE KINEMATICS PER REGION 

Most of the extinction regions contain multiple dense cores, thus the relationship be- 
tween the ensemble of dense cores and the extinction region can also be examined. In 
particular, the core-to-core (vlsr) motions can be compared to the spread of motions seen 
over the region as a whole in an "environmental" tracer (e.g., 13 CO); this should reveal how 
connected the dense cores are to the lower density gas in the region. If the dense cores are 
connected to the large-scale non-thermal motions, then the dispersion in the core centroid 
velocities would be expected to be similar to the global velocity dispersion of the region, as 
measured in 13 CO. If, however, the core-to-core motions are much smaller than the global 
velocity dispersion of the region, then this may be an indication that the dense cores have be- 
come detached from the large-scale motions within the cloud, arising from turbulent flows, 
for example. Alternatively, the spatial distribution of cores may lead to them tracing a 
smaller volume of material than the 13 CO. 
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5.1. Method of Calculation 

Section 4.1 describes the calculation of the extinction region velocity dispersions. There 
are several approaches that can be taken to quantify the core-to-core motions within each 
extinction region. The first, method A, is a purely statistical measure, which involves simply 
taking the standard deviation of the centroid velocities of all of the dense cores within each 
region, av met h.A- This is effective when the spectrum of each core is well described by a 
single, relatively narrow, velocity component. In cases where these conditions do not hold, 
interpretation of the standard deviation can be more difficult. This method also does not 
account for individual core linewidths, which could represent a significant fraction of the total 
dense gas velocity dispersion of the core gas. The second method, method B, accounts for 
both the core-to-core motion and the internal velocity dispersion through summing all of the 
spectra and measuring the resultant total velocity dispersion, av met h.B- This method may 
be biassed if some cores are much brighter than others, and can be particularly challenging 
for molecules with hyperfi ne trans itions, due to different components blending together in 



the summed spectrum. In IKJT07I . we adopted method A for measuring the N 2 H + relative 
motions and method B for measuring the C 18 relative motions; here, we again adopt method 
A for N2H + , but show that method B yields similar results. Note that any velocity gradient 
across the region is implicitly included in both measures and is not accounted for separately; 
the same is true for the extinction region veloctiy dispersion. 



5.2. Results 

Figures [7] through [12] qualitatively show our main result (examined in more detail be- 
low) - the dense cores have a much lower core-to-core velocity dispersion than the velocity 
dispersion of the 13 CO regions they reside in. In the top panel of each figure, a comparison 
is shown between the cumulative 13 CO spectrum across the entire region with measures of 
the dispersion seen in all the N2H + dense cores. The dashed light green line shows the result 
for N2H + from method A above (the standard deviation of the N2H + centroid velocities, 
centered on the mean centroid velocity), while the dark green line shows the result from 
method B above (the dispersion of the sum of the N 2 H + spectra). As well, the dark blue 
line shows the result from method B applied to the C 18 spectra in the region. All of these 
cumulative spectra are compared to the region as a whole in 13 CO (black and grey shading). 
The figures show that the core-to-core velocity dispersion is similar when measured using 
either method A or method B, and that this velocity dispersion is significantly smaller than 
the velocity dispersion of the region as a whole. 



We now determine the above result more rigorously. In 



KJT07 



, we preferred the use 



-12- 



of method A to measure the N2H + core-to-core velocity dispersion, due to the difficulties 
associated with method B for N 2 H + outlined in Section 5.1. In that work, only starless 
cores were considered, due to the potential for protostars to have inherited more motion 
from surrounding sources. Here, we measure the standard deviation both for the starless 
cores only as well as the full sample of cores (the full sample was plotted in Figures [7] to 
H2l) . Table [3] provides the values derived for the core-to-core velocity dispersion for both 
method A (both starless cores and the full sample) and for method B (full sample only). 
Figure d3] shows a comparison of the core-to-core velocity dispersion measured by method 
A for the starless cores ((Jv met h.A,si) and all cores (crv me th.A,aii) versus the velocity dispersion 
of the region (cri 3 co,reg)- As can be seen in the figure, the core-to-core velocity dispersion 
tends to be similar in both samples; variations between the two may be attributed at least in 
part to small number statistics. The magnitude of variation expected due to small number 
sampling is further discussed in the following subsection. 

The standard deviation of core centroid velocities measured in each extinction region 
is a good description of the width of the distribution, despite the small number of cores 
typically present. In the regions with a sufficient number of cores, we tested this in two 
ways. First, we ran KS tests and found that the distribution of velocities in each region 
was consistent with being drawn from a normal distribution with a median and standard 
deviation equal to the value measured, with over 80% probability in most cases. As a second 
test, we compared the standard deviation to the median absolute deviat ion, a robust me asure 



of the width of distributions, i.e., little-effected by outlying points (e.g.. lAndersenll2008l ). For 
a normal distribution, the median absolute deviation is 1.48 times smaller than the standard 
deviation; using this conversion factor, we find the standard deviation as inferred from the 
median absolute deviation agrees to within 12%. Based on these tests, it is reasonable to 
assume that the standard deviation is a good descriptor of the width of the core centroid 
velocity dispersion within each region. 

As Figure [13] suggests, there is no apparent variation with the core-to-core velocity dis- 
persion with the region velocity dispersion. In fact, all are consistent with being originating 
from the same parent distribution. Two-sided KS tests of each pair of regions yields proba- 
bilities of 65% in all but one case (and 12% in the remaining case) that the core velocities 
are drawn from the same sampl^f). 



5 All regions were compared with region 7, since that region has the largest number of cores. The median 
core velocity in each region was subtracted from all the core velocities prior to comparison to account for 
the large-scale velocity gradient across Perseus. 
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5.2.1. Effect of Sampling a Small Number of Cores 

The difference in size between the core-to-core and 13 CO regional velocity dispersion 
is sufficiently large that it cannot be attributed to the small number of cores within each 
region except for the regions with a very small number of cores. We ran a set of simulations 
to determine the size of the error induced by using a small number of cores to measure the 
underlying total velocity dispersion. Scaling the regional velocity dispersion to a value of 
one, we selected a sample of N cores objects with velocities randomly drawn from a normal 
distribution (using IDL's randomn function), and calculated the standard deviation of the 
object-to-object velocity. We ran 10000 simulations for each value of N cores . Figure E] shows 
the results of this calculation. The squares indicate the mean scaled velocity dispersion de- 
termined for the 10000 simulations for every value of N cores , while the vertical lines indicate 
the 68.2% confidence levels, corresponding to the 1 a level in a normal distribution. Sam- 
pling a very small number of cores tends to increase the difference between the measured and 
intrinsic velocity dispersion, but does not systematically lower the measured velocity disper- 
sion by a significant amount when more tha n a few cores a re present, in agreement with the 



full analytic result from statistics (see, e.g., iKenneyl 119391 ); note that the full expression for 



the bias in the standard deviation is more complex than the oft-quoted yjn/{n— 1)0 and is 
important to use for very low number samples. 

The diamonds in Figure [Til show the observed core-to-core velocity dispersions (including 
both the starless cores and protostars), scaled by the 13 CO velocity dispersion of the region. 
This figure clearly shows that all of the observed core-to-core velocity dispersions lie well 
below the range of values expected for a population of objects sharing the same velocity 
dispersion as the parent 13 CO gas. Weighting the ratios between the core-to-core and region 
velocity dispersions after correction for the bias in the standard deviation gives a value of 
0.5 ± 0.1, for either the full sample, or when restricted to the regions with good spectral 
coverage. This implies that if all regions have similar ratios of core-to-core versus regional 
velocity dispersions, a ratio of around one half is a reasonable estimate. 

We conclude from this analysis that the dense cores tend to have a much smaller core- 
to-core motion than is present in the mean 13 CO gas, roughly half the value, and that this 
result is not due to small number statistics. 



12 r(%) 

3 The full correction factor is <\ / r ~ , „ — T~ 
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5.2.2. Effect of Spatial Sampling 

Since the dense cores tend to be preferentially found concentrated within a small por- 
tion of the extinction region, it is reasonable to ask whether the difference in size of the 
core-to-core and total region velocity dispersions is caused by spatial sampling rather than 
a physical difference between the two. For example, the dense cores could be at the centre 
of the potential well of the region and thus naturally be expected to have a smaller veloc- 
ity dispersion than the mean gas spanning the entire potential well. Here, we show that 
the core-to-core velocity dispersions are intrinsically smaller than the typical 13 CO velocity 
dispersions, even accounting for the difference in area sampled on the plane of the sky. 

The bottom panel of Figures [7J through [12] show the total 13 CO spectrum for each 
extinction region with differing areas considered. The black line and grey shading (identical 
to the one in the bottom panel) represents the sum of all 13 CO spectra within the region. 
The purple line shows the result that would be obtained if only the spectra within the area 
spanned by the cores were considered. In order to obtain this estimate, we summed only 
the spectra that lie within the minimal box (in RA and dec) containing all of the dense 
cores within the region. The resultant spectrum is similar to the total (black) spectrum, and 
certainly is not similar in width to the core-to-core result (green lines on the top panels). 

We also considered an even more restrictive case - the red line in the bottom panel of 
Figures [71 through [121 shows the resultant spectrum if the 13 CO spectra only at the locations 
of the dense cores are used in the summation. This again shows marked similarity to the 
spectrum for the region as a whole, although with a lower signal to noise level. The 13 CO 
spectra for each region when restricted to the area spanned by or locations of the dense 
cores can be somewhat narrower than the 13 CO spectrum over the entire region, however, 
this difference is not sufficient to explain the smaller-still distribution of core-to-core velocity 
dispersions within each region. We measured the ratio in the width (Gaussian sigma) fit to 
the areally-restricted 13 CO spectra versus the spectra for the entire region, and find values 
of 88 ± 11% and 86 ± 16% for the spectra for the region spanned by cores and at only the 
core locations, respectively. Since the core-to-core velocity dispersion is approximately 50% 
of the total region velocity dispersion, the bulk of the difference cannot be explained by the 
clustered locations of dense cores. 



5.3. Interpretation 

Our analysis indicates that the core-to-core velocity dispersion is significantly smaller 
than the velocity dispersion of the lower density 13 CO gas in the larger extinction regions 
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in which the cores are found. Although this result is robust against the potential biasses 
that we have investigated, interpretation is hampered by a lack of knowledge of the full 
3D structure of the cloud. If, for example, the 13 CO within each extinction region can be 
thought of as belonging to a single, large, coherent entity, then the fact that the core-to-core 
velocity dispersions are significantly smaller than that of the entire region would indicate 
that the dense cores are more bound than the region as a whole. On the other hand, if the 
13 CO gas actually comprises a number of smaller, independent structures along and across 
the LOS, then the 13 CO may have a larger velocity dispersion than the cores because of the 
multiple structures probed along the LOS. Even if the 13 CO gas comprises a single entity 
along the line of sight, if the dense cores probe only a fraction of this mate rial, a smalle r 
velocity dispersion would be expected. Larson's linewidth-size relationship ( ILarsonlll98ll ). 
for example, would predict a velocity dispersion half as large for the cores as the region if the 
cores traced a length scale roughly one quarter the length traced by 13 CO. We do not have 
any way to determine the three dimensional cloud structure, hence the ratio in core-to-core 
vs region velocity dispersion we measure cannot have a a definitive interpretation. It is, 
however, a straightforward benchmark that can be used to evaluate star formation models. 



6. REGIONAL VELOCITY GRADIENTS 

Finally, we analyze the large-scale features in the velocity structure of each extinction 
region. We calculate the gradient of 13 CO centroid velocities across each extinction region 
using the 5' convolve d 13 CO data (see Secti on 4.1) and applying the gradient calculation 



method formulated in iGoodman et al.l (119931 ) . Figures [151 and [161 show the centroid velocity 
measured for each cell in the convolved 13 CO map (colourscale) and the gradient measured 
(white arrows). The overall regional velocity gradient is shown in the bottom left corner of 
each plot. Table H] provides the magnitude and orientation of the gradient measured in each 
region, along with the formal error to the gradients fit; note that these errors are often small 
due to the large number of points fit. 



6.1. Cores Relative to the Gradient 

If the dense cores are connected to the large-scale motions within each region, then 
they would be expected to follow the same large-scale gradient as is observed in 13 CO. We 
determine the velocity expected at each position in the region by using the velocity at the 
centre of the extinction region as a zero-point and extrapolating outwards using the gradient 
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measured, i.e., 

Vfiow,mod(x, y) = v 13CO (xc, y c ) + G -d (2) 

or 

Vflow,mod(x,y) = V 13C o(Xc,y C ) + | Q | &\\ (3) 

where visco is the velocity observed at the central position (x c ,y c ), Q is the gradient 
(in km s -1 pc -1 ), d is the distance from (x,y) to (x c ,y c ), and d\\ is the distance from (x,y) 
to (x c ,y c ) in the direction parallel to the gradient's direction. 

Figure [T71 shows the difference in core centroid velocity, Vnw+, and the velocity inferred 
at the core's position from the regional velocity gradient, Vfi ow ^ mo d as a function of RA. As in 
Figure [2], the background greyscale image shows the PV diagram for 13 CO, with the value at 
each RA shifted so the peak intensity is at zero. A similar scatter in the velocity difference 
is found in all extinction regions. 

The left panel of Figure [17] shows the distribution in velocity difference between the 
cores and the regional flow, VN2H+ — ^siow^mod for the entire cloud (red), as well as the 
velocity difference between individual 13 CO pointings and the regional flow, Vi 3 co — Vfiow,mod 
(black). While the cores do not perfectly follow the gradient-determined velocity, they tend 
to have a deviation similar to that of the typical 13 CO gas; the standard deviations of both 
VN2H+ - v f i OWt mod and v 13C o - Vfi OWtmod are 0.45 km s _1 . 



6.2. Large-scale Energetics 

The magnitude of the gradient across the region can be compared to the total velocity 
dispersion measured for the region. Regions dominated by large-scale motions such as shear 
or rotation, would be expected to have a high ratio of the velocity dispersion inferred from the 
gradient to the total velocity dispersion, whereas regions dominated by small-scale random 
motions would be expected to have a lower ratio. 

We approximate the velocity dispersion inferred from the gradient across each region 
as the gradient multiplied by the diameter of the region. Since the 13 CO emission across 
each region is, to first order, constant, this is a (somewhat) better estimate than using, for 
example, the radius of the region; all locations along the region contribute approximately 
equally to the total region spectrum. (Note that this estimate would not be appropriate 
for centrally concentrated objects.) Figure [T8l shows a comparison of the velocity dispersion 
inferred from the gradient across each region versus the total dispersion observed in 13 CO. 
The ratio of large-scale versus total velocity dispersion is around 50% for the regions with 
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well-determined measurements. A few of the regions show ratios above 100%, which can 
be understood by the fact that our estimate of the large-scale velocity dispersion is crude. 
For example, a varying column density across a region would tend to weight a limited range 
in velocities, and hence decrease the total velocity dispersion measured, while leaving the 
large-scale velocity dispersion determined from the gradient unchanged. 

Turbulent motions that are domina t ed by large-scale modes can also produce an observ- 
able gradient. iBurkert fc Bodenheimerl (120001 hereafter BB00) investigated the magnitude 
of the gradi ent expected from large -scale turbulent modes for cores similar t o the NIF> cores 
observed in 
BG98), and 



Goodman et al 



Goodman et al. 



1993, hereafter G93), iBarranco fc Goodman! ( 119981 . hereafter 
( 119981 ). BB00 ran several thousand realizations of a 3D turbu- 
lent velocity field with power spectra of P(k) oc k~ A and P(k) cx k~ 3 , where P(k) is the power 
at wavenumber k\ the oft-adopted Kolmogorov (incompressible turbulence) power spectrum 
falls within this range with P(k) oc k^ 11 ^ 3 . Along with this turbulent velocity field, BB00 
adopted a spherically symmetric, centrally-condensed density distribution. They then cal- 
culated the resultant spectra observab le in 2D, and determined the gradient across the core 
using the same method as described in lG93l (and also used in our analysis). The P(k) oc k~ 4 
turbulent power spectrum has more power in the largest modes, and, as expected, was shown 
to generally produce larger observed gradients. 

Although the BB00 models were run under the assumption of a single size scale (cores 
with a radius of 0.1 pc), BB00 as sume a Larson-type scaling law with a oc -R ' 5 to compare 
their results to that of the entire lG93l sample of NH 3 cores. Figure [TH shows the velocity 
dispersions and diameters of the extinction regions we measure as well as the G93 NH3 
cores. The extinction regions follow a similar trend to the G93 cores, hence the BB00 
models can be compared to the extinction regions in the same manner as was done for the 
G93 cores. In Figure [2DJ we show the BB00 predictions for gradients resulting from turb ulent 
power spectra with P(k) oc k~ 4 and P(k) oc k~ 3 (diagonal lines), as well as the lG93l data 
(triangles) and our own 13 CO measurements (diamonds). For comparison, we also show a 
recent measurement of the size and yelocity gradient o f an embedded prot ostar (IRS 1) within 
the isolated Bok globule, BHR 71 (jChen et all 120081 ). To plot the lG93l d ata, we adopt the 
same procedure as BB00 - the diameter used is the geometric mean of the lG93l FWHM sizes 
of the major and minor axes, and updated gradient measurements for the seven sources in 
BG98I replace those from lG93l . The solid diagonal lines indicate the most likely values of the 
gradient found in the BB00 simulations, while the dotted lines indicate the half maxi mum 
likelihood values. While the extinction regions we analyze are larger than the cores in lG93l . 
our observations clearly follow the same trend predicted by the BB00 model. Although 
the centrally concentrated density distribution adopted in the BB00 model is not a good 
representation of the density distribution in the regions we analyze, BB00 state that their 
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results are not strongly dependent on the density distribution - they found similar results 
when a constant density distribution was adopted. The extension of their predictions to the 
13 CO gradients in the extinction regions measured here is therefore reasonable. A generalized 
formulatio n of the velocity gradient f or various density profiles and turbulent power spectra 



Kratter fc Matznerl (120061 ). and also shows similar results are expected regardless 



is given m 
of the assumed density distribution. 



7. DISCUSSION - ENERGETICS 

The above analyses all point towards a common picture of the energetics o f dense cores 



within the context of the molecular clouds they inhabit. Previous work (e.g.. IKJT07I ) has 
shown that dense cores have little motion with respect to their immediate surrounding gas. 
Our analysis comparing the N2H + and 13 CO data has shown that the small motions continue 
to larger scales. Within the larger extinction regions that the cores inhabit, the cores have 
smaller core-to-core velocity dispersions than the typical 13 CO gas. Interpretation of this 
observation is hampered by the lack of knowledge of the full 3D structure of the lower density 
cloud material. The lower density gas traced by 13 CO may be tracing multiple structures 
(or a larger length scale) along the line of sight that the cores are not associated with. 

The large-scale low density gas appears to be dominated by turbulent motions. Across 
each extinction region, the gradients in 13 CO centroid velocity were measured and appear to 
account for roughly half of the total velocity dispersion of the region. This is approximately 
the same size as is predicted to arise from large-scale modes of a turbulent power spectrum 
with P(k) oc k~ 4 (BB00). The velocities of the dense cores within each extinction region 
follow the velocity gradient found in the 13 CO gas. This is consistent with a picture where 
both the dense cores and 13 CO gas obey Larson's Laws; the cores trace a smaller volume of 
gas than the 13 CO both across the plane of the sky and along the line of sight, and so show a 
smaller core-to-core velocity dispersion than the regional 13 CO velocity dispersion, while still 
tracing the same large-scale velocity gradient. A direct comparison of the ratio in the core- 
to-core and extinction region velocity dispersions with the prediction from Larson's Laws 
could be made once similar measurements exist for several molecular clouds. Determining 
this effect quantitatively for a single region is complicated by the lack of knowledge of the 
line of sight distribution of cores and 13 CO. With a large enough sample of regions, however, 
this effect should be quantifiable across the ensemble. 

Comparison of the data to expectations from virial equilibrium is unclear at best. As 
usual, the surface terms of the virial equation are unknown and must be excluded. In 
addition, the gas within each extinction region must be assumed to be a single coherent 
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entity both across the sky and along the line of sight and the substantial velocity gradient 
across each region discussed above must also be ignored. With these many assumptions, the 
virial velocity dispersion of each region could then be estimated in some of the traditional 
ways - using the mass and size of each region (as in KTJ07) or using the 13 CO velocity 
dispersion as a proxy (similar to Larson 1981). The core-to-core motion within each region 
would then be interpreted as "sub- virial" , since their velocity dispersion is much smaller than 
either of these two estimates of the region's virial velocity dispersion. Given the questionable 
assumptions that are required to make the virial argument, in particular ignoring the large- 
scale gradient across each region and the possibility of multiple structures along the line of 
sight, we do not recommend interpreting the data under this framework. Instead, future 
'observations' of numerical simulations will be able to reveal the range of physical conditions 
that can match our results. 



8. CONCLUSIONS 



We analyze the kinematic relationship between dense cores an d the larger environments 
they inha bit using pointed N2H" 1 " and C 18 Q ob servations (IKJT07I ) as well as a 13 CO spec- 
tral cube ( IRidge et al.ll2006t iPineda et al.ll2008l ) of the Pe rseus m olecular cloud. We present 
a series of kinematic measures complementing those in (1KJT071 ) that can be used to test 
models of star formation. We choose to investigate physically-motivated, but procedurally 
straightforward measures of kinematics in Perseus. Through this, we hope that future simu- 
lations and as well as observations of other star-forming regions will adopt similar measures 
to allow constraints of models of star formation across differing environments. 

We find the dense cores have small motions with respect to the 13 CO gas along the 
same line of sight, with a standard deviation of 0.33 km s" 1 , or 0.39 times the 13 CO velocity 
dispersion at that location, while the sound speed in the ambient medium is ~ 0.23 km s _1 . 
The small motions we measure are an upper limit on the true motion between the N 2 H + 
and 13 CO, since the 13 CO data have a beamsize roughly twice the size of the N 2 H + data. 
The 13 CO may also be tracing additional kinematic structures that are not traced by N2H + . 
The core-to-core motions within the larger extinction regions they inhabit are about half 
of the total ( 13 CO) velocity dispersion within the region. The dense cores tend to follow 
the overall velocity gradient across each extinction region to the same extent as the 13 CO 
gas. The gradients across each extinction region tend to account for roughly half of the 
total measured velocity dispersion, and are of the magnit ude as expected by large-scal e 
turbulent modes following a power spectrum of P(k) cx k~ A (IBurkert fc Bodenheimerll2000l ). 
The kinematic measures presented here should be relatively straighforward to obtain both 
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in other large observational datasets as well as numerical simulations of star formation. 
Together with the measures presented in IKJT07I . we provide a set of benchmarks by which 
future star formation models can be evaluated. 
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A. GLOSSARY OF KINEMATIC MEASURES 

A large number of kinematic measures are presented in this paper in order to provide 
a variety of diagnostics on the molecular cloud. For clarity, Table |5] summarizes all of these 
kinematic measures, along with a brief description of their meaning. 
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Table I. N 2 H+ to CO Relative Motions 



C ls O a 13 CO point b 13 CO region d 



1N2-H- 
# 


VN2H+ - VC180 
(km s ) 


(km s J 


VN2H+ — V13CO 
(km s J 


o"13CO 
(km s ) 


CjOOQ 

Fit? 


Region 
# 


VN2H+ — Vl3CO,reg 

(km s J 


<?13CO,reg 

(km s J 


Fit? 


4 


0.13 


0.91 


0.11 


0.58 


Y 


1 


0.19 


0.64 


Y 


5 


0.01 


1.00 


0.21 


0.64 


Y 


1 


0.15 


0.64 


Y 


6 


0.19 


0.86 


0.35 


0.67 


Y 


1 


0.32 


0.64 


Y 


8 


-0.13 


0.53 


-0.19 


0.50 


Y 


1 


0.02 


0.64 


Y 


12 


-0.10 


0.40 


0.32 


0.69 


Y 


1 


0.71 


0.64 


Y 


13 


0.01 


0.41 


0.44 


0.67 


Y 


1 


0.94 


0.64 


Y 


15 


-0.17 


0.86 


-0.09 


0.61 


Y 


2 


0.14 


0.91 


Y 


16 


0.06 


1.36 


0.50 


0.78 


Y 


2 


1.07 


0.91 


Y 


18 


0.33 


1.23 


0.54 


0.83 


Y 


2 


1.12 


0.91 


Y 


19 


-0.01 


0.67 


0.03 


0.64 


Y 


2 


-0.25 


0.91 


Y 


20 


0.02 


0.62 


0.14 


0.65 


Y 


2 


-0.23 


0.91 


Y 


21 


0.20 


1.00 


0.41 


0.71 


Y 


2 


0.22 


0.91 


Y 


22 


-0.49 


1.13 


-0.55 


0.68 


Y 


2 


-0.69 


0.91 


Y 


23 


0.02 


0.75 


-0.19 


0.82 


Y 


2 


-0.53 


0.91 


Y 


24 


-0.02 


0.85 


0.17 


0.75 


Y 


2 


-0.06 


0.91 


Y 


25 


0.12 


0.81 


0.32 


0.64 


Y 


2 


0.26 


0.91 


Y 


26 


0.04 


0.57 


-0.22 


0.67 


Y 


2 


-0.24 


0.91 


Y 


27 


-0.48 


0.57 


-0.54 


0.55 


Y 


2 


-0.56 


0.91 


Y 


27 


-0.02 


0.57 


-0.08 


0.55 


Y 


2 


-0.10 


0.91 


Y 


28 


-0.21 


0.98 


-0.42 


0.57 


Y 


2 


-0.34 


0.91 


Y 


30 


0.10 


0.44 


0.19 


0.56 


Y 


2 


0.04 


0.91 


Y 


33 


0.11 


0.50 


0.65 


0.79 


N 


2 


0.65 


0.91 


Y 


34 


-0.03 


0.58 


-0.20 


0.63 


Y 


2 


-0.59 


0.91 


Y 


36 


0.03 


0.35 


0.26 


0.61 


Y 


2 


0.59 


0.91 


Y 


41 


0.06 


0.39 


0.39 


0.89 


Y 


2 


-0.35 


0.91 


Y 


53 


0.10 


0.50 


0.36 


0.65 


Y 


3 


0.73 


1.01 


Y 


55 


0.06 


0.50 


0.57 


0.90 




3 


0.45 


1.01 


Y 


63 


-0.06 


0.64 


-0.08 


0.70 


Y 


4 


-0.35 


0.88 


Y 


67 


-0.40 


1.00 


-0.66 


0.71 


Y 


5 


-0.46 


0.89 


Y 


68 


0.11 


1.04 


0.31 


0.92 


Y 


5 


0.23 


0.89 


Y 


71 


0.25 


0.85 


-0.07 


0.85 


Y 


5 


-0.03 


0.89 


Y 


72 


-0.20 


1.20 


-0.38 


0.94 


Y 


5 


-0.42 


0.89 


Y 


73 


-0.22 


0.95 


-0.24 


0.81 


Y 


5 


-0.26 


0.89 


Y 


74 


-0.09 


0.86 


-0.18 


0.70 


Y 


5 


0.15 


0.89 


Y 


75 


-0.07 


0.89 


-0.16 


0.83 


Y 


5 


-0.06 


0.89 


Y 


76 


-0.12 


0.73 


-0.20 


0.61 


Y 


5 


-0.18 


0.89 


Y 


77 


-0.04 


0.79 


-0.15 


0.84 


Y 


5 


-0.05 


0.89 


Y 


78 


-0.01 


1.09 


-0.19 


0.91 


Y 


5 


-0.02 


0.89 


Y 


79 


0.25 


1.02 


0.09 


0.91 


Y 


5 


0.13 


0.89 


Y 


84 


-0.21 


0.79 


0.23 


0.96 


Y 


5 


-0.08 


0.89 


Y 


85 


-0.09 


0.98 


-0.16 


0.71 


Y 


5 


-0.28 


0.89 


Y 


86 


0.01 


0.96 


0.25 


0.74 


Y 


5 


0.21 


0.89 


Y 


87 


0.00 


0.41 


-0.48 


0.85 


Y 


5 


-0.79 


0.89 


Y 


90 


-0.07 


1.08 


-0.19 


0.94 


Y 


6 


0.41 


0.91 


Y 



Table 1 — Continued 







C 1S a 


13 CO point b 






13 CO 


region d 




N 2 H+ a 


VN2H+ — 


v ClSO CT C180 


VN2H+ - V13CO °"13CO 


Good 


Rcgion c 


V N2H+ — v 13CO,reg 


u lZCO,reg 


Good 


# 


(km s 


- 1 ) (km s- 1 ) 


(km s _1 ) (km s _1 ) 


Fit? 


# 


(km s _1 ) 


(km s _1 ) 


Fit? 



91 


0.05 


0.56 


0.68 


0.87 


Y 


- 


- 


- 


N 


92 


0.02 


0.57 


- 


- 


- 


6 


-0.44 


0.91 


Y 


93 


-0.23 


0.60 


_ 


- 


- 


6 


-0.58 


0.91 


Y 


94 


-0.08 


0.67 


- 


- 


- 


6 


-0.62 


0.91 


Y 


95 


0.05 


0.52 


-0.07 


0.39 


Y 


7 


0.53 


1.25 


Y 


96 


0.01 


0.71 


-0.10 


1.17 


N 


7 


-0.08 


1.25 


Y 


97 


-0.01 


0.49 


-0.00 


0.48 


Y 


7 


-0.12 


1.25 


Y 


98 


0.08 


1.05 


-0.29 


0.91 


Y 


7 


-0.15 


1.25 


Y 


99 


0.01 


0.39 


-0.09 


0.74 


N 


7 


-0.08 


1.25 


Y 


99 


0.00 


0.46 


-0.72 


0.74 


N 


7 


-0.71 


1.25 


Y 


100 


-0.03 


0.71 


-0.14 


0.81 


Y 


7 


-0.10 


1.25 


Y 


101 


-0.14 


1.42 


0.27 


1.15 


Y 


7 


0.59 


1.25 


Y 


102 


0.12 


1.11 


0.18 


1.04 


N 


7 


0.22 


1.25 


Y 


103 


0.19 


0.86 


-0.58 


1.22 


Y 


7 


-0.83 


1.25 


Y 


103 


0.55 


0.86 


-0.22 


1.22 


Y 


7 


-0.47 


1.25 


Y 


104 


0.05 


1.41 


0.88 


1.57 


Y 


7 


0.97 


1.25 


Y 


106 


0.19 


1.92 


0.06 


1.19 


Y 


7 


-0.13 


1.25 


Y 


107 


0.25 


0.41 


-0.38 


1.22 


Y 


7 


-0.38 


1.25 


Y 


109 


-0.04 


1.49 


0.62 


1.51 


Y 


7 


0.86 


1.25 


Y 


110 


0.41 


1.88 


0.43 


1.26 


Y 


7 


0.35 


1.25 


Y 


111 


0.14 


1.59 


-0.52 


1.25 


Y 


7 


-0.80 


1.25 


Y 


111 


0.35 


2.13 


0.40 


1.25 


Y 


7 


0.12 


1.25 


Y 


112 


0.22 


1.51 


0.65 


1.38 


Y 


7 


0.78 


1.25 


Y 


113 


0.24 


1.46 


0.36 


1.18 


Y 


7 


0.35 


1.25 


Y 


115 


-0.30 


1.53 


-0.20 


1.22 


Y 


7 


-0.46 


1.25 


Y 


116 


-0.17 


1.53 


0.04 


1.01 


Y 


7 


0.21 


1.25 


Y 


118 


-0.28 


1.31 


-0.27 


1.21 


Y 


7 


-0.28 


1.25 


Y 


118 


0.12 


1.31 


0.13 


1.21 


Y 


7 


0.12 


1.25 


Y 


121 


-0.06 


0.92 


0.20 


0.94 


Y 


7 


-0.40 


1.25 


Y 


122 


-0.08 


0.96 


-0.29 


0.74 


N 


7 


0.36 


1.25 


Y 


123 


-0.15 


0.94 


-0.18 


0.72 


N 


7 


0.57 


1.25 


Y 


124 


0.13 


1.57 


0.06 


0.98 


Y 


7 


-0.58 


1.25 


Y 


125 


-0.07 


0.70 


-0.46 


0.93 


Y 


7 


-0.34 


1.25 


Y 


126 


-0.02 


1.64 


-0.17 


1.12 


Y 


7 


-0.78 


1.25 


Y 


127 


0.04 


0.43 


-0.02 


0.46 


Y 


9 


0.18 


0.87 


Y 


128 


0.07 


1.40 


-0.28 


1.13 


N 


7 


-0.41 


1.25 


Y 


129 


0.01 


0.23 








9 


-0.80 


0.87 


Y 


130 


-0.01 


0.67 








9 


-0.51 


0.87 


Y 


131 


0.05 


0.38 








9 


-0.58 


0.87 


Y 


132 


-0.37 


1.01 








9 


-0.59 


0.87 


Y 


133 


-0.16 


1.07 








9 


-0.30 


0.87 


Y 


134 


-0.04 


0.70 








9 


-0.49 


0.87 


Y 


135 


0.04 


0.70 








9 


-0.64 


0.87 


Y 


136 


0.24 


1.34 








9 


-0.72 


0.87 


Y 
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Table 1 — Continued 





C 18 a 






13 CO point b 






13 co 


region d 




N 2 H+ a 


"JV2H+ — VC1SO 


0"C180 


VN2H+ — VriCO CT 13CO 


Good 


Region 


VN2H+ ~ '"13CO,reg 


fl3CO,reg 


Good 


# 


(km s" 1 ) 


(km s" 1 ) 


(km s 1 


) (km s" 1 ) 


Fit? 


# 


(km s" 1 ) 


(km s — 1 ) 


Fit? 


136 


0.04 


0.45 








9 


0.58 


0.87 


Y 


138 


-0.17 


0.90 








9 


-0.29 


0.87 


Y 


139 


0.17 


0.97 








9 


-0.27 


0.87 


Y 


140 


0.00 


0.27 








9 


-0.02 


0.87 


Y 


143 


0.00 


0.46 


0.12 


0.56 


Y 


9 


-0.27 


0.87 


Y 


145 


-0.18 


1.31 


-0.42 


0.60 


N 


10 


-0.94 


0.81 


Y 


146 


0.15 


0.88 


0.06 


0.71 


Y 


10 


0.13 


0.81 


Y 


147 


0.09 


0.86 


0.05 


0.70 


Y 


10 


0.28 


0.81 


Y 


148 


-0.25 


0.53 


0.35 


0.80 


Y 


10 


0.75 


0.81 


Y 


148 


0.22 


1.02 


0.27 


0.80 


Y 


10 


0.67 


0.81 


Y 


149 


0.16 


1.07 


0.01 


0.79 


N 


10 


0.15 


0.81 


Y 


150 


-0.04 


0.68 


-0.20 


0.74 


N 


10 


-0.31 


0.81 


Y 


151 


0.11 


0.56 


0.08 


0.47 


Y 


11 


0.23 


0.88 


Y 


152 


0.03 


0.74 


-0.04 


0.62 


N 


10 


-0.28 


0.81 


Y 


155 


0.04 


0.49 


-0.20 


0.45 


N 


11 


-0.14 


0.88 


Y 



a N2H+ core number, N2H+ to C ls O relative velocity, and C ls O velocity dispersion from lKirk. Johnstone fc Tafallal i 20071) . 

b N2H" 1 " to 13 CO relative velocity, 13 CO velocity dispersion at the point of the N2H + observation, and whether the 13 CO fit was deemed 
good. Only good fits were used in the analysis. Dashes denote regions where 13 CO data does not exist. 

c Associated extinction region. 

d N2H+ to 13 CO relative velocity, 13 CO linewidth for the region and whether the 13 CO fit was deemed good. Only good fits were used 
in the analysis. 
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Table 2. Region Velocity Dispersions in 13 CO 



Region a 


Areal a 


rr b 

u reg, Gauss 


Good b 


&reg,FWQM ° 


# 


Cov. (%) 


(km s _1 ) 


Fit? 


(km s _1 ) 


1 


100 


0.64 


Y 


0.60 


2 


100 


0.91 


Y 


0.88 


3 


95 


1.01 


N 


0.98 


4 


55 


0.88 


N 


1.06 


5 


100 


0.89 


Y 


0.86 


6 


52 


0.91 


Y 


0.86 


7 


100 


1.25 


N 


1.22 


8 


100 


0.90 


Y 


0.84 


9 


40 


0.87 


Y 


0.82 


10 


100 


0.81 


Y 


0.74 


11 


84 


0.88 


Y 


0.84 



a Extinction region number and fractional coverage of the area observed in 13 CO over the 
extinction region. 

b The Gaussian-fit 13 CO region velocity dispersion and whether the fit was judged to be 
good. 

c The equivalent 13 CO region velocity dispersion as measured from the FWQM. See text 
for details. 
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Table 3. Core-to-Core Velocity Dispersions 



Region 


crv sh A a 


<7V Mi a a 




# 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


1 


0.39 


0.36 


0.35 


2 


0.55 


0.54 


0.31 


3 


0.20 


0.20 


0.20 


4 






0.14 


5 


0.28 


0.28 


0.21 


6 


0.09 


0.49 


0.20 


7 


0.54 


0.50 


0.56 


8 
9 


0.32 


0.38 


0.41 


10 


0.55 


0.56 


0.35 


11 


0.26 


0.26 


0.26 



a Core-to-core velocity dispersion derived 
for starless cores and all cores using method 
A. Dashes denote regions with less than two 
cores. See text for details. 

b Core-to-core velocity dispersion derived 
for all cores using method B. Dashes denote 
regions with no cores. See text for details. 
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Table 4. 13 C0 Gradients Across Each Extinction Region 



Region a Q b Error b Angle c 


Error c 


D d 


1 G 1 D e 


# (km s _1 pc _1 ) (km s _1 pc _1 ) (degrees) 


(degrees) 


(pc) 


(km s _1 ) 



1 


0.1456 


0.0009 


84.3 


0.3 


1.88 


0.27 


2 


0.1115 


0.0002 


-118.78 


0.08 


2.71 


0.30 


3 


0.1845 


0.0007 


84.5 


0.2 


1.79 


0.33 


4 


0.811 


0.002 


-132.5 


0.1 


1.63 


1.32 


5 


0.2252 


0.0004 


-82.98 


0.07 


1.40 


0.32 


6 


0.667 


0.005 


23.4 


0.3 


1.10 


0.73 


7 


0.3230 


0.0003 


-6.11 


0.04 


2.16 


0.70 


8 


0.827 


0.003 


-48.6 


0.1 


1.10 


0.91 


9 


0.934 


0.005 


-56.1 


0.2 


1.09 


1.01 


10 


0.420 


0.002 


-42.3 


0.2 


0.94 


0.39 


11 


0.627 


0.002 


-33.4 


0.1 


0.75 


0.47 



a The extinction region number 

b The velocity gradient and associated error measured across each extinction region. See 
Section 6 for more detail. 

c Orientation of the gradient and associated error. The angle is given in degrees clockwise 
from North. 



d Extinction region diameter from IKJD06 



s The velocity dispersion inferred from the gradient across each region. 



Table 5. Summary of Kinematic Measures Used 



Symbol Formula Description 



v x Centroid velocity of species x 

a x Velocity dispersion (Gaussian sigma) of species x 
VN2H+ — v x Relative motion between N2H + and species x 

Cnorm,x VN2H a + x ~ Vx Normalized relative motion 

c r eg,Gauss Region velocity dispersion using a Gaussian fit a 

&reg,FWQM FW QM reg / Region velocity dispersion measured using the full width at quarter maximum a 

<Jv me th.A a ( v ) Dispersion of core centroid velocities within an extinction region ("method A") 

<7v me th.B Velocity dispersion of the cores within an extinction region ("method B") b 

D reg Diameter of extinction region 

Q Velocity gradient across a region 

( x c,Uc) position of centre of region used in gradient calculation 

d\\ Distance between (x, y) and (x c , y c ) parallel to the gradient 

Vfiow,mod(x,y) v(x c ,y c )+ I Q | d|| Velocity inferred from the bulk gradient 

&reg,mod G x D reg Region velocity dispersion inferred from the bulk gradient 



a Measured using the cumulative spectrum of the entire extinction region 
b Measured using the cumulative core spectra within the extinction region 



I I I I I I I I I I I I I I I I I I I I I I I I I I I 




Fig. 1. — The greyscale image shows the 13 CO integrated intensity map from the FCRAO 
telescope. The overlaid black contours show extinction measurements derived from 2MASS 
data at levels of Ay = 3, 5, and 7 mag. The red circles show the locations of all of the 
IRAM 30 m N 2 H + observations. The extinction regions used in the analysis here are denoted 
by the coloured contours. See text for details. 



31 




Fig. 2. — The distribution of centroid velocity difference between N 2 H + , C 18 0, and 13 C0 
as a function of RA (left panel), and the full distribution (right panel). In the left panel, 
the blue diamonds indicate the distribution of VN2H+ — vciso, while the red circles indicate 
the distribution of VN2H+ — ^i3CO- The empty diamonds indicate locations where there is 
no good Gaussian fit to the 13 CO data (or 13 CO observations do not exist), while the filled 
diamonds indicate where there is a good Gaussian fit to the 13 CO data. The greyscale 
image shows the position- velocity diagram for the 13 CO gas and the grey contours show the 
half-maximum value at every RA. At each RA, the cumulative spectrum has been shifed so 
that the peak intensity falls at a value of zero on the diagram. The greyscale ranges from 
2 K km s _1 (black) to K km s" 1 (white). The approximate RA of well-known star forming 
regions are indicated by the text. In the right panel, the dotted blue line indicates the full 
distribution of v^2H+ ~ v ciso, while the solid blue line shows only the points where there is 
a good Gaussian fit to the 13 CO data. The distribution of VN2H+ — ^i3CO is shown in red. 
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Fig. 3. — The distribution of normalized velocity differences, (norm, between N 2 H + , C O, 
and 13 CO as a function of RA (left panel). The total distribution is shown in the right panel. 
See Figure [2] for details on the plotting conventions used. 
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Fig. 4. — A comparison of the different measures of the velocity dispersion within each ex- 
tinction region, illustrating that the measures are similar in regions where the 13 CO coverage 
is good. The 13 CO velocity dispersion is measured with both a single Gaussian (diamonds) 
and the equivalent Gaussian sigma derived from measuring the full width at quarter max- 
imum (FWQM) of the spectra (squares). Red points denote the single Gaussian fits that 
were judged by eye to fit well, while blue points denote regions where the 13 CO data covers 
80% or more of the extinction region, i.e., the spectrum is a good representation of that 
for the entire region. In all regions where there is good 13 CO coverage (blue points), both 
velocity dispersion measurements (diamonds and squares) are similar, regardless of whether 
the single Gaussian fits were judged to be good (red points) or not. The remainder of this 
paper uses the single Gaussian fit measures, restricted to regions where there is good 13 CO 
coverage. 
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Fig. 5. — The distribution of the centroid velocity difference between N 2 H + and the 13 CO 
region in which the core is located as a function of RA (left panel), and the total distribution 
(right panel). The greyscale image in the left panel shows the position- velocity (PV) diagram 
for the 13 CO gas with the grey contours showing the half-maximum value at every RA. At 
each RA, the cumulative spectrum has been shifed so that the peak intensity falls at a value 
of zero on the diagram. The greyscale ranges from 2 K km s _1 (black) to K km s _1 (white). 
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Fig. 6. — The distribution of normalized centroid velocity differences, ( n0 rm between the 
N2H + and the region in which the core is located as a function of RA (left panel) and the 
total distribution (right panel). See previous figure for the plotting conventions used. 
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Fig. 7. — A comparison of cumulative spectra within extinction region 1. Top panel: The 
13 CO spectrum of the entire region from the top panel (black line and grey shading), the 
C 18 spectrum summed over all locations in the region with C 18 detections ("method B", 
blue line), the resultant fit for a similar summation for N 2 H + ("method B", dark green line), 
and the standard deviation of the centroid N 2 H + velocities in each region ("method A"), 
centered on the mean value (light green dashed line). Bottom panel: 13 CO spectra derived 
from summing the spectra over the entire region (black line), summing over the box in which 
the N 2 H + dense cores are found only (purple dot-dashed line), and summing over only the 
points where where N 2 H + dense cores were detected (red line). Note the 13 CO spectral data 
has a full range of ~ —10 to 30 km s _1 ; the truncated range of to 12 km s" 1 plotted covers 
the range of velocities where N 2 H + and C 18 emission was observed. 
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Fig. 8. — Cumulative spectra within extinction regions 2 and 3. See Figure [TJfor the plotting 
conventions used. 
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Fig. 9. — Cumulative spectra within extinction regions 4 and 5. See Figure [TJfor the plotting 
conventions used. Note that in extinction region 4, only one dense core was observed, so 
the standard deviation of N 2 H + centroid velocities (dashed green line) cannot be plotted in 
the top panel. The purple dashed and red solid line are identical in the bottom plot for the 
same reason. 
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Fig. 10. — Cumulative spectra within extinction regions 6 and 7. See Figure [7J for the 
plotting conventions used. The offset in N2H + centroids in extinction region 6 is caused by 
the small number of cores (four) in the region, one of which is faint enough that it impacts 
the measurement using method A only. 




Fig. 11. — Cumulative spectra within extinction regions 8 and 9. See Figure [7J for the 
plotting conventions used. Note that there were no N2H + detections in extinction region 8. 
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Fig. 12. — Cumulative spectra within extinction regions 10 and 11. See Figure [7] for the 
plotting conventions used. 



-40- 




0.4 



0.6 0.8 1.0 

Region Velocity Dispersion (km s"') 



1.2 



1.4 



Fig. 13. — The spread in core-to-core velocity dispersions found in each extinction region 
versus the velocity dispersion of that region. The diamonds show the values derived using 
all of the dense cores within each region, while the circles show the values derived using only 
the starless cores within each region. The open symbols denote the values derived for the 
extinction regions with poor 13 CO coverage. The dotted line shows a 1-1 relationship. 
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Fig. 14. — The expected range of observed core-to-core velocity dispersions relative to the 
distribution they were drawn from for a given number of cores. The boxes and vertical 
lines show the scaled mean and standard deviation of the core-to-core velocity dispersion for 
varying numbers of cores observed. The diamonds show our observed values scaled to the 
velocity dispersion of the region, with the open symbols denoting the extinction regions with 
poor 13 CO coverage. The dotted horizontal line indicates an observed velocity dispersion 
equal to that of the parent sample, while the dashed line shows the scaled mean core-to-core 
velocity dispersion values multiplied by one half. 
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Fig. 15. — The overall velocity gradient measured for each extinction region in 13 CO for 
extinction regions 1-6. The colourscale indicates the centroid velocity of each cell, while the 
arrows indicate the direction of the local gradient. The overall gradient is indicated in the 
bottom left corner of each panel. Diamonds show the positions of the N 2 H + cores in each 
region. 
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Fig. 16. — The overall velocity gradient measured for each extinction region in 13 CO for 
extinction regions 7 to 11. The colourscale indicates the centroid velocity of each cell, while 
the arrows indicate the direction of the local gradient. The overall gradient is indicated in 
the bottom left corner of each panel. Diamonds show the positions of the N 2 H + cores in 
each region. 
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Fig. 17. — The distribution of centroid velocity differences between the dense cores and 
the velocity inferred at their location from the 13 CO velocity gradient measured in their 
extinction region (points) as a function of RA (left panel). The RAs of well-known star- 
forming regions are noted, with the exception of L1455 where the 13 CO observations do not 
cover a large enough area to allow for a model of the gradient to be applied. As in Figure [2} 
the background greyscale image shows the PV diagram of the 13 CO gas and the grey contours 
show the half-maximum value at every RA, with the spectrum at each RA shifted so that 
the peak emission falls at 0. The total distrbution of VN2H+ — Vfiow, mo d is shown in the right 
panel in red, and for comparison, the total distribution of visco — Vfi ow ^ mo d is shown in black. 
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Fig. 18. — The ratio of large-scale to total velocity dispersion observed in each extinction 
region. The open symbols show the results for the three extinction regions where the 13 CO 
coverage is less than 80% and hence less reliable. 
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Fig. 19. — The observed velo city dispersion versus diameter for all of the extinction regions 
(diamonds) and NH 3 cores in lGoodman et al.l (119931 ) (triangles). The open symbols show the 
results for the regions where the 13 CO coverage was less than 80%. The dotted line shows 
the slope of the scaling relationship assumed by iBurkert fc Bodenheimerl (120001). a oc R 5 . 
Note that a scale factor for the relationship is not given in lBurkert fc Bodenheimerl ( 120001 ). 
so the line plotted is a guide to the slope only. 
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Fig. 20. — The gradient measured for each extinction region versus the diameter of the region 
(diamonds) . The open diamonds show the results for the three extinction regions where the 
13 CO coverage was less than 80% and h ence less reliable. The triangles show the da t a from 
the NH3 cores of iGoodman et al.l (119931 ). while the cross shows data from I Chen et al.l ( 120081 ); 
see text for details. The red diagonal lines show the range (dotted) and most likely values 
(solid) predicted by BBOO for a k~ 4 turbulent power spectrum, while the blue lines show the 
range (dotted) and most likely values (solid) for a k~ 3 turbulent power spectrum. 



